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We investigate the evolution, final fate, and nucleosynthetic yields of rotating and 
non-rotating very massive stars (VMS) of zero metallicity. First we address the issue of 
mass loss during hydrogen burning due to vibrational instabilities. We find that these 
objects are much more stable than what was found in previous studies of VMS of solar 
composition, and expect only negligible mass loss driven by the pulsations. As these 
^ ■ stars thus reach the end of their evolution with massive helium cores, they encounter the 
pair-creation instability. We find that for helium core masses of ~ 64 . . . 133 M Q these 
stars are completely disrupted with explosion energies of up to ~ 10 53 erg and eject up 
to ~ 60 M Q of 56 Ni. Stars with more massive helium cores collapse into black holes. We 
present the first calculations that follow the collapse of such a massive rotating star and 
predict that X-ray burst and significant gravitational wave emission could result. 

... 

^ . 1. INTRODUCTION 

Recently, three-dimensional cosmological simulations have reached sufficient resolution 
on small scales to begin to address the star formation problem of the first generation of 
so-called Population III (Pop III) stars Though the formation and nature of the 

first generation of metal-free stars have been investigated since over thirty years 0-0, and 
many others] , the extent to which they differed from present day stars in ways other than 
composition is still debated. While considerable uncertainty remains, these simulations 
suggest that the first generation of stars may have been quite massive, ~ 100 — 1000 M 
giving us motivation to examine the properties, evolution, and fate of such stars. 



2. VIBRATIONAL (IN-) STABILITY AND MASS LOSS 

Above a critical mass, main sequence stars are vibrationally unstable due to the destabi- 
lizing effect of nuclear reactions in their central regions (e- mechanism) |]TT],|T2|. According 
to non-linear calculations, such an instability leads to mass loss rather than catastrophic 



disruption fi3-15|. Previous investigations of stars of solar-like compositions or slightly 



metal-poor stars [|T|,|l^-|2(| indicated that stars above a few 100 M would lose a substan- 



tial amount of mass during hydrogen burning. However, the structure of metal-free stars 
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Figure 1. Post-explosive nucleosynthetic yields of helium cores as a function of mass 
(black lines; left axis). 20 Ne and 36 Ar have yields similar to those of respectively 24 Mg 
and 40 Ca. The thick gray line gives the explosion energy in "foe" (1 foe = 10 51 erg, about 
the explosion energy of a "typical" core collapse supernova; right axis) . 



is significantly different from those that even contain a trace of initial metals. We have 
thus performed, to our knowledge, the first stability analysis of zero metallicity stars with 
100 M < M < 500 M Q . 

We find that the different structure and the higher central temperatures of these stars 
make them less vibrationally unstable, i.e., the timescale for the growth of the amplitude 



is much longer, and thus, based on the non-linear calculations of [13|, we derive consid- 



erably lower mass loss rates than for more metal rich counterparts, and find that even 
our 500 M star loses only a negligible amount of mass during central hydrogen burning. 
During central helium burning only stars with initial masses above ~ 300 M become pul- 
sationally unstable during the last few 10,000 yr of their evolution. Therefore we estimate 
that stars of < 500 M can even retain their hydrogen envelope until they encounter the 
pair creation instability. A more detailed description of our analysis can be found in [El 



3. NUCLEOSYNTHETIC YIELDS OF PAIR CREATION SUPERNOVAE 



The evolution of very massive stars has already been followed by many authors before 
[^, e.g.]. We followed the evolution and nucleosynthesis of rotating and non-rotating VMS 
Pop III stars from onset of hydrogen burning through the electron-positron pair creation 
instability until one year after the explosion f23[, using recent stellar "input physics" and 
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extended nuclear reaction networks. 

In Fig. [I] we give the results of calculations of non-rotating plain helium cores evolved 
without mass loss. It displays the yields of the ejecta as a function of helium core mass 
and the resulting explosion energy. For a 12 C(a,7) 16 rate of 1.7x that of |23|, in the 
extreme case of a 133.3 M helium core we find an explosion energy of about 9.4xl0 52 erg 
and ejection of 57 M Q of 56 Ni, almost half of all ejecta. The ejecta from these objects 
are dominated by "a-nuclei", but some nickel isotopes are made in the a-rich freeze- 
out in the center of the most massive cores. Essentially no elements above the iron 
group are produced. Below the iron group, the elements of even charge number are 
produced in about solar abundance ratios, while elements of odd mass number are strongly 
underproduced. 

Below ~ 64 M we do not find prompt explosions of the helium cores, but after a strong 
pair-instability induced pulse the outer layers are ejected while the central parts of the 
star fall back and contract again and are either disrupted in subsequent pulses or evolve 
through to iron core collapse pi| and could become a collapsar, similar to the models of 
||26|| . Above ~ 133 M the photo-disintegration of heavy elements into a-particles, and 
a-particles into nucleons, consumes so much energy that the center of the star continues 
to collapse into a black hole (see below). Only in the mass range shown in Fig. [1|, which 
corresponds to initial stellar mass of ~ 150 . . . 250 M Q , the star is completely disrupted in 
one pulse. 

4. SIMULATIONS OF THE COLLAPSE TO A BLACK HOLE 

We have followed the collapse of a rotating 300 M stars that had a rotation of 10 % 
Keplerian on the zero-age main sequence and which formed a ~ 180 M helium core. The 
ID model was mapped into a 2D Lagrangian SPH code when the central density exceeded 
5 x 10 10 g cm -3 and the further evolution was followed including, for the first time, neutrino 
trapping during the collapse. Sufficient angular momentum was present for triaxial de- 
formations to grow significantly in the "proto-black hole" , with the possible consequence 
of a strong gravity wave signal. After the formation of the initial black hole, it accreted 
mass to ~ 140 M at a rate of 10 . . . 100 M s _1 . At this point, a centrifugally supported 
disk forms that may produce jets along the poles through magneto-hydrodynamical ef- 
fects [J27| . Though the efficiencies of the mechanisms under consideration are still very 
speculative, such a jet could lead either to an energetic, jet-driven explosion or an x-ray 
transient, depending on whether the hydrogen-rich envelope of the star is still present or 
has been lost to a possible binary companion star. A more detailed description is given 



5. SUMMARY AND CONCLUSIONS 

Recent results on the formation of Pop III stars indicate that very massive stars of zero 
metallicity could have formed. Our analysis of the pulsation properties indicates that stars 
of < 500 M are far less vibrationally unstable than previously though and thus may reach 
the end of their evolution as objects with massive helium cores that encounter the pair- 
creation instability. This has interesting consequences for the nucleosynthetic yield from 
these objects, and may even allow for the observation of their explosions despite the high 
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red shift. Above an initial mass of ~ 300 M Q the core collapses into a black hole and may 
cause an X-ray burst or a strong gravity wave signal that may be accessible to experiments 
in the near future. 

The survival of these objects as very massive stars until the end of their evolution also 
has interesting consequences for the re-ionization of the early universe, as these VMS have 
an effective surface temperature of ~ 10 5 K [ B5 , ^T |, which means that they very efficiently 



turn nuclear energy into ionizing photons. 
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